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ABSTRACT 

Context. Blue Stragglers Stars (BSSs) are thought to form in globular clusters by two main formation channels: i) mergers induced 
by stellar collisions and if) coalescence or mass-transfer between companions in binary systems. The detailed study of the BSS 
properties is therefore crucial for understanding the binary evolution mechanisms, and the complex interplay between dynamics and 
stellar evolution in dense stellar systems. 

Aims. We present the first comparison between the BSS specific frequency and the binary fraction in the core of a sample of Galactic 
globular clusters, with the aim of investigating the relative efficiency of the two proposed formation mechanisms. 
Methods. We derived the frequency of BSSs in the core of thirteen low-density Galactic globular clusters by using deep ACS@HST 
observations and investigated its correlation with the binary fraction and various other cluster parameters. 

Results. We observed a correlation between the BSS specific frequency and the binary fraction. The significance of the correlation 
increases by including a further dependence on the cluster central velocity dispersion. 

Conclusions. We conclude that the unperturbed evolution of primordial binaries could be the dominant BSS formation process, at 
least in low-density environments. 
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1. Introduction 

Blue Straggler Stars (BSSs) are objects that, in the color- 
magnitude diagram (CMD) of evolved stellar populations, lie 
along an extension of the Main Sequence (MS), in a region 
which is brighter and bluer than the Turn-off (TO). First dis- 
covered by Sandage (1953) in M3, they have been observed in 
all Galactic globular clusters (GCs; Piotto et al. 2004), in the 
field population (Carney et al. 2005), and in dwarf galaxies of 
the local group (Momany et al. 2007). 

Their location in the CMD suggests that BSSs have masses 
of 1.2 4- 1.5 M Q , significantly larger than those of normal stars in 
old stellar systems (like GCs). Thus, they are thought to have in- 
creased their mass during their evolution. Two mechanisms have 
been proposed for their formation: ;) the merger of two stars in- 
duced by stellar collision (COL-BSSs; Hills & Day 1976) and 
ii) coalescence or mass-transfer between two companions in a 
binary system (MT-BSSs; McCrea 1964). The two formation 
channels are thought to act with different efficiencies accord- 
ing to the cluster structural parameters (Fusi Pecci et al. 1992) 
and they can work simultaneously within the same cluster in dif- 
ferent radial regions, corresponding to widely different stellar 
densities (Ferraro et al. 1997; Mapelli et al. 2006). Indeed, col- 
lisions are more frequent in the central region of GCs, because 
of the high stellar density, while MT-BSSs mainly populate the 
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cluster periphery, where binary systems can more easily evolve 
in isolation without suffering exchange or ionization due to grav- 
itational encounters. The whole scenario is further complicated 
by the cluster dynamical evolution that leads massive systems 
(like binaries and BSSs) to sink toward the cluster center in a 
time-scale comparable to the cluster relaxation time. 

A possible tool for distinguishing COL-BSSs from MT- 
BSSs is based on high-resolution spectroscopic analysis. In fact, 
anomalous chemical abundances are expected at the surface of 
BSSs resulting from mass-transfer activity (Sarna & de Greve 
1996), while they are not predicted for COL-BSSs (Lombardi 
et al. 1995). However, such studies have just become feasible 
and they are limited to only a small number of BSSs in just one 
cluster (47 Tucanae; Ferraro et al. 2006). 

As an alternative way for getting insights on the relative ef- 
ficiency of the two formation mechanisms, here we investigate 
possible correlations between the BSS population and the host 
cluster properties. 

2. Observations and Data Reduction 

The photometric database considered here is the same used in 
Sollima et al. (2007; hereafter S07) for measuring the frequency 
of binary systems in the core of thirteen low-density Galactic 
GCs. It consists of a set of high-resolution images, obtained with 
the ACS on board HST, through the F606W (V 606 ) and F814W 
ih\A) filters. As described in S07, the target clusters are char- 
acterized by high Galactic latitude (b > 15°), low reddening 
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(E(B - V) < 0.1) and low projected central density (logpg < 5, 
where p' Q is in units of M arcrnin -2 ). These selection crite- 
ria have been chosen in order to limit the effects of field con- 
tamination, differential reddening and crowding on the determi- 
nation of the binary fractions. The selected target clusters and 
their main physical parameters are listed in Table [T] The central 
luminosity density po and the absolute V band magnitude My 
are from Djorgovski (1993), the age fq is from Salaris & Weiss 
(2002 fH and the projected central velocity dispersion cr v is from 
McLaughlin & Van der Marel (2005). A detailed description of 
the data reduction and calibration procedure can be found in S07 
together with the resulting CMDs for all the target clusters. 

As discussed in S07, the binary sequence is well defined and 
distinguishable from the MS and the fraction of binary systems 
(ijbin) found within the cluster core radius ranges from about 10 
to 50% (the value of are listed in Table[T]l. Moreover, a num- 
ber of BSSs populating the bright part of the CMD is also evi- 
dent in the data sample. In the following section we describe the 
adopted procedure to derive the BSS frequencies in the core of 
the target clusters. 



3. BSS population selection 

Our primary criterion for the definition of the BSS sample is 
based on the location of stars in the (/gi4, V(,o6 - hu) CMD. 
Since the considered GCs have different distances, reddening, 
ages and metallicities, in order to define a homogeneous criterion 
for the BSS selection, we shifted in magnitude and color the 
CMD of each cluster to match that of M55 (see also Ferraro 
et al. 1995). For minimizing the contamination from sub-giant 
branch, horizontal branch (HB) and field stars, we limited the 
analysis to colors bluer than ( V606 _ hi4) < 0-47 and magnitudes 
in the range 15.1 < 7si4 < 18.1. The adopted BSS selection box 
(the same for all clusters) is shown in Figure[T]and the number of 
selected BSSs (A^bss) in eacn target cluster is reported in Table 
ffl 

In order to compare the frequency of BSSs in different GCs, 
the size of the total cluster population must be taken into con- 
sideration. Commonly, the number of BSSs is normalized to 
the number of stars belonging to a given evolutive sequence 
(usually the HB, or the Red Giant Branch, RGB; Ferraro et al. 
1995; Piotto et al. 2004; Leigh et al. 2007). Unfortunately, in 
some clusters of our sample the number of HB and RGB stars 
is very small, which would make the derived BSS frequency 
highly uncertain. We therefore normalized the number of se- 
lected BSSs, to the number of MS stars (Nms) in the magni- 
tude range 18.1 < Igu < 19.1 and with colors comprised within 
A(V606 ~ ^8m) < 0.05 from the cluster MS mean ridge line. Stars 
in this magnitude range are bright enough to have the same com- 
pleteness level of BSSs (<p > 95%) and cover a similar mass 
range in all clusters: following the theoretical stellar isochrones 
by Cariulo et al. (2004) and taking into account the differences 
in age and metallicities among the selected GCs, the TO mass 
varies between 0.8 M Q < M < 0.9 M . Hence, we considered 
Nms a good indicator of the size of the sampled cluster popula- 
tion. 

The possible effect of background and foreground field star 
contamination on the BSS and MS selections has been estimated 
by using the Galaxy model of Robin et al. (2003). For this pur- 
pose a synthetic catalog covering an area of 0.5 square degree 



Table 1. Main physical parameters of the target globular clusters 



Name 


log Po 


M v 


t 9 






Nbss 


F 








Gyr 


Kms 1 


% 




% 


NGC 288 


1.80 


-6.63 


11.3 


2.79 


11.6 


43 


1.9 


NGC 4590 


2.52 


-7.73 


11.2 


3.59 


14.2 


62 


1.5 


NGC 5053 


0.51 


-7.07 


10.8 


1.79 


11.0 


15 


1.0 


NGC 5466 


0.68 


-6.83 


12.2 


1.86 


9.5 


39 


1.7 


NGC 5897 


1.32 


-7.27 


12.3 


2.94 


13.2 


53 


1.3 


NGC 6101 


1.57 


-6.82 


10.7 


2.91 


15.6 


60 


1.1 


NGC 6362 


2.23 


-6.69 


11.0 


3.83 


11.8 


42 


1.5 


NGC 6723 


2.71 


-7.95 


11.6 


5.52 


16.1 


35 


0.5 


NGC 6981 


2.26 


-7.16 


9.5 


3.66 


28.1 


65 


1.5 


M55 


2.12 


-7.47 


12.3 


3.73 


9.6 


32 


0.8 


Arp2 


-0.35 


-5.30 


9.7 


0.72 


32.9 


33 


2.3 


Terzan 7 


1.97 


-5.99 


7.4 


2.55 


50.9 


33 


3.4 


Palomar 12 


0.68 


-4.78 


6.4 


0.76 


40.8 


9 


2.5 



centered on each cluster (with the centers taken from Djorgovski 
& Meylan 1993) has been retrieved, both in the V and the / 
Johnson-Cousin bands. A sub-sample of stars scaled to the size 
of the ACS field of view (202" x 202") has been randomly 
extracted, and the Johnson-Cousin magnitudes have been con- 
verted into the ACS photometric system by using the transfor- 
mations of Sirianni et al. (2005). The number of contaminating 
field stars turns out to be neglegible (<2) in all the clusters of our 
sample. Then, we counted the number of field stars included in 
the BSS and MS selection boxes (N^ e f d and N^ u , respectively), 
and we computed the BSS specific frequency as: 



(1) 



The resulting values for the thirteen target clusters are listed 
in Table [T] In order to check the suitability of the adopted nor- 
malization criterion, we compared the derived BSS specific fre- 
quencies with those measured by Leigh et al. (2007, and upt- 
dated by Sills 2008, private communication) for the four clusters 
in common with our sample (namely, NGC 4590, NGC 6362, 
NGC 6723 and NGC 6981). These authors normalized the num- 
ber of BSSs in the cluster cores to that of red giant branch stars 
in the same magnitude range. In spite of the different normal- 
ization criterion, the specific frequencies calculated for the four 
clusters in common correlate very well (with a correlation coef- 
ficient r = 0.92). 

4. Correlations between BSS specific frequency 
and cluster parameters 

With the aim of shading light on the BSS formation mechanisms, 
we correlated the BSS specific frequency (F ) calculated above 
with the core binary fraction (£bm) obtained by S07 and with 
various other physical parameters of the thirteen GCs of the sam- 
ple. In particular, we searched for possible correlations with the 
cluster central density, total magnitude, age, and central velocity 
dispersion (see Table [TJ, as well as with the stellar collision rate 
(defined as F = p 2 , r^/cr v ; Pooley & Hut 2006), the global metal- 
licity ([M/H], from Ferraro et al. 1999^}, the cluster evaporation 
rate (v, from Gnedin & Ostriker 1997), concentration (c, from 



1 Since NGC 6981 is not included in the list of Salaris & Weiss 
(2002), we converted the ages measured by De Angeli et al. (2005) into 
the Salaris & Weiss (2002) scale. 



2 For NGC 6101, NGC 6362, NGC 6723 and Palomar 12 not in- 
cluded in the list of Ferraro et al. (1999) we transformed the metallicity 
[Fe/H] from Zinn & West (1984) into the global metallicity [M/H] fol- 
lowing the prescriptions of Ferraro et al. (1999). 
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Fig.l. (7 8 i4, y 606 - hu) CMD of M55. The adopted selection 
boxes for the BSS and MS populations are shown. The fiducial 
BSSs are marked with open circles. 



Trager et al. 1995), half-mass relaxation time (f r /,) and half-mass 
radius (r/„ from Djorgovski 1993). 

The Pearson's linear correlation coefficients and confidence 
levels for the considered parameter pairs are listed in Table [2] 
The strongest correlation is found between the BSS specific fre- 
quency and the binary fraction. A noticeable correlation with the 
absolute magnitude and anticorrelations with the cluster age and 
central velocity dispersion are also apparent. Considering that 
the age estimates are quite uncertain and span a narrow range 
of values in our sample, any firm conclusion about the possible 
dependence of F on the cluster age is premature at the moment. 
However, if confirmed, such an anticorrelation might be a conse- 
quence of the analogous anticorrelation between binary fraction 
and age already discussed by S07, and would suggest that binary 
disruption processes in the core of GCs become more efficient 
with time, thus reducing the fraction of both binaries and BSS 
in the core of the oldest clusters (as suggested by Ivanova et al. 
2005, but at odds with Hurley et al. 2007). In any case, even if 
the anticorrelation between F (or £/,,„) and tg is not confirmed, 
the correlation between the BSS specific frequency and the bi- 
nary fraction still holds. The latter is shown in the upper panel 
of Figure [2] and it suggests that the main formation channel of 
BSS in these clusters is the unperturbed evolution of primordial 
binary systems. 

We stress that, since the considered parameters do not con- 
stitute an orthonormal basis, most of the less significant cor- 
relations could be artificially induced by primary correlations 
among other parameters. In order to investigate the possible de- 
pendence of F on more than one cluster parameter, we applied 
the Bayesian Information Criterion (BIC) test (Schwarz 1978) 
to our dataset. We assumed the BSS specific frequency F to be a 
linear combination of a sub-sample of p parameters (Aj) selected 
among those listed above: 



7 fU 



-2> 



and we calculated the quantity 
P 



BIC = log L n 



log N, 



(3) 



where L p is the likelihood and N is the dimension of the sample 
(N = 13 in our case). The likelihood is calculated as log L p = 
Zf=i log P p j, where 



p,j 



cr F 



V2lr 



(4) 



Ai + a 



(2) 



with Fj being the BSS fraction of cluster j, and cr F being the 
residual of the fit. The p parameters that maximize the quantity 
BIC are the most probable correlators with F. For our dataset 
the maximum value of BIC is obtained with p = 2, A\ = t;un 
and A2 — o~ v . The correlation between F, and cr v is shown in 
Figure|2] together with the best fit line. All the other 3-variate or 
higher-order correlations turn out to be non-significant. 

A note of caution is worth: as can be seen, these correla- 
tions are essentially due to the four clusters (namely Terzan 7, 
Palomar 12, Arp 2 and NGC 6981) with the largest (> 20%) bi- 
nary fraction, while the other GCs in our sample define a group 
with an average binary content ~ 10% and a BSS frequency 
F ~ 1.2%. In particular, all the correlations that we have found 
are mainly driven by Terzan 7, Palomar 12 and Arp 2, which are 
the GCs not only with the largest binary and BSS fractions, but 
also with the lowest absolute magnitudes and the lowest ages. 
These three clusters are also the most distant from the Sun and 
they are thought to belong to the Sagittarius Stream (Bellazzini 
et. al 2003). Thus, they might be stellar systems with intrinsi- 
cally different origins and properties, and the correlations found 
here might not be appropriate for the overall class of "genuine" 
Galactic GCs (GGCs). By excluding these three clusters from 
the sample, the anticorrelation with the total luminosity still 
holds (even if with lower significance), in agreement with what 
already found by Piotto et al. (2004), Sandquist (2005) and Leigh 
et al. (2007) for larger samples of GCs with different charac- 
teristics, and for BSS populations selected and normalized with 
different criteria. Instead, the anticorrelation with the age disap- 
pears, thus bringing our results in agreement with those of Leigh 
et al. (2007; see also De Marchi et al. 2006). The correlation 
between the BSS and the binary fractions also disappears ex- 
cluding the three "Sagittarius" clusters (SgrGCs), with the range 
in spanned by the remaining sample being very small. Note 
however that both the average BSS and binary frequencies are 
higher in SgrGCs «F BS s) = 2.7% and <£,,■„> = 41.5%) than 
in GGCs «F B ss> = 1-3% and = 14%). Thus, even as- 

suming that SgrGCs were born with different initial conditions 
from GGCs, the observational facts presented here allow to draw 
(at least) the "conservative" conclusion that binary-rich envi- 
roments tend to produce more BSS in the low-density cluster 
regime (where collisions are expected to play a minor role in the 
formation/destruction of binaries). 

5. Conclusions 

We measured the BSS specific frequency in the core of thirteen 
low-density Galactic GCs and investigated its correlation with 
different dynamical and general cluster parameters. We found 
evidences that, at least in this density regime, binary-rich envi- 
ronments are more efficient in producing BSS. No correlations 
have been found with the cluster central density, concentration, 
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Table 2. Pearson's correlation coefficients (r) and confidence 
levels (P) 



Parameter 




P 


£bin 


0.820 


0.9998 


My 


0.774 


0.9985 


tg 


-0.758 


0.998 


cr v 


-0.594 


0.98 


[Af/fl] 


0.403 


0.92 


Po 


-0.311 


0.85 




0.299 


0.84 


r 


-0.158 


< 0.75 


c 


0.129 


< 0.75 


V 


0.099 


< 0.75 


log frfi 


0.020 


< 0.75 




binaries is to increase their binding energy (i.e. decrease their 
orbital separation; Heggie 1975), this implies that low velocity 
dispersion GCs should host a larger fraction of hard (and close) 
binaries, able to both survive possible stellar encounters, and ac- 
tivate mass-transfer and/or merging processes between the com- 
panions. A larger fraction of BSSs formed by the evolution of 
primordial binaries is therefore expected in lower velocity dis- 
persion GCs (see also Davies et al. 2004). Such an effect of <x v 
(in terms of both hardening and shrinking the binary systems) 
might be, in turn, at the origin of the inverse correlation between 
the BSS frequency and the cluster total luminosity (mass) ob- 
served in open clusters (De Marchi et al. 2006), low density 
GCs (Sandquist 2005), as well as high density GCs (Piotto et 
al. 2004; Leigh et al. 2007). Indeed, the more massive GCs have 
larger central velocity dispersions (as a consequence of the virial 
theorem for systems with similar radii, as GCs; see Fig. 1 of 
Djorgovski 1995). For most of these clusters, however, the bi- 
nary fraction is still unknown and the trivariate correlation be- 
tween F, ^bi„ and <x v cannot be derived. If its significance and 
its interpretation in terms of the velocity dispersion effect is con- 
firmed also in high-density clusters, then stellar collisions might 
play a secondary role in the production of BSSs, and the evolu- 
tion of primordial binaries should always be the dominant pro- 
cess. For the moment, however, these conclusions remain spec- 
ulative, since the sample of GCs analyzed here, in spite of be- 
ing the largest to date with known binary fraction, is still too 
small for statistically reliable assessments. Enlarging the sample 
of GCs with known binary and BSS fraction is therefore essen- 
tial and urgent to verify the findings presented in this paper on a 
more robust statistical basis. 
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Fig. 2. Specific frequency of BSSs as a function of the core bi- 
nary fraction (upper panel) and of the best-fit linear combi- 
nation of core binary fraction and cluster velocity dispersion 
(bottom panel, obtained with a = -4.62). Solid lines indicate 
the best-fit lines; dashed lines mark the boundaries of the +0.005 
range in F with respect to the best-fit lines. 



stellar collision rate, and half-mass relaxation time, in agreement 
with the results of Piotto et al. (2004) and Leigh et al. (2007). 
These evidences indicate that the collisional channel for the BSS 
formation has a very small efficiency in low-density GCs, while 
the mechanisms involving the unperturbed evolution of binary 
systems are dominant. 

The higher significance of the trivariate correlation among 
the BSS frequency, the binary fractions and the cluster velocity 
dispersion, indicates that, for a given binary fraction, the BSS 
specific frequency decreases with increasing velocity dispersion. 
This finding might be connected with the effect of the cluster ve- 
locity dispersion in the dynamical evolution of binary systems. 
In fact, a small cluster velocity dispersion corresponds to a lower 
energy limit between soft and hard binaries^] i.e., to a larger frac- 
tion of hard binary systems. Since the natural evolution of hard 
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